home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
IRIX 6.2 Development Libraries
/
SGI IRIX 6.2 Development Libraries.iso
/
dist
/
complib.idb
/
usr
/
share
/
catman
/
p_man
/
cat3
/
complib
/
SPOSL.z
/
SPOSL
Wrap
Text File
|
1996-03-14
|
2KB
|
67 lines
SSSSPPPPOOOOSSSSLLLL((((3333FFFF)))) SSSSPPPPOOOOSSSSLLLL((((3333FFFF))))
NNNNAAAAMMMMEEEE
SPOSL - SPOSL solves the real symmetric positive definite system A * X
= B using the factors computed by SPOCO or SPOFA.
SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
SUBROUTINE SPOSL(A,LDA,N,B)
DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
On Entry
AAAA REAL(LDA, N)
the output from SPOCO or SPOFA.
LLLLDDDDAAAA INTEGER
the leading dimension of the array A .
NNNN INTEGER
the order of the matrix A .
BBBB REAL(N)
the right hand side vector. On Return
BBBB the solution vector X . Error Condition
AAAA division by zero will occur if the input factor contains a zero on the
diagonal. Technically, this indicates singularity, but it is usually
caused by improper subroutine arguments. It will not occur if the
subroutines are called correctly and INFO .EQ. 0 . To compute
INVERSE(A) * C where C is a matrix with P columns
CALL SPOCO(A,LDA,N,RCOND,Z,INFO)
IF (RCOND is too small .OR. INFO .NE. 0) GO TO ...
DO 10 J = 1, P
CALL SPOSL(A,LDA,N,C(1,J)) 10 CONTINUE LINPACK. This version dated
08/14/78 . Cleve Moler, University of New Mexico, Argonne National Lab.
Subroutines and Functions BLAS SAXPY,SDOT
PPPPaaaaggggeeee 1111